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Abstract 

Recent work from authors across disciplines has made substantial con- 
tributions to counting rules (Maxwell type theorems) which predict when 
an infinite periodic structure would be rigid or flexible while preserving 
the periodic pattern, as an engineering type framework, or equivalently as 
an idealized molecular framework. Other work has shown that for finite 
frameworks, introducing symmetry modifies the previous general counts, 
and under some circumstances this symmetrized Maxwell type count can 
predict added finite flexibility in the structure. 

In this paper we combine these approaches to present new Maxwell type 
counts for the columns and rows of a modified orbit matrix for structures 
that have both a periodic structure and additional symmetry within the 
periodic cells. In a number of cases, this count for the combined group 
of symmetry operations demonstrates there is added finite flexibility in 
what would have been rigid when realized without the symmetry. Given 
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that many crystal structures have these added symmetries, and that their 
flexibility may be key to their physical and chemical properties, we present 
a summary of the results as a way to generate further developments of both 
a practical and theoretic interest. 

Key Words: framework rigidity, periodic, symmetry, crystal systems, orbits 



1 Introduction 



The theory of periodic frameworks has undergone rapid and extensive develop- 
ment in the last four years ( |Borcea an d Streinu 2010a|b Malestein and Theran 



2010} Ross 2010} . We now have necessary conditions (call them Maxwel 



counts) for such frameworks to be rigid, either with a fixed lattice of transla- 
tions or with a flexible lattice of translations. Underlying much of the recent 
work are finite 'lattice rigidity matrices' for the equivalence classes of vertices 
and edges under the infinite group of translations Z d in d-space. With the 
corresponding count of periodicity-preserving trivial motions under these con- 
straints (typically d translations), the number of rows, e, and columns, dv + I 
(where / is the number of lattice parameters) of these 'orbit matrices' lead to nec- 
essary Maxwell type counts for a framework to be infinitesimally rigid (Borcea 
and Streinu] |2010bipalestein and Theran|[2^T0||Ross|[201u) : e>dv + l-d. 

The theory of finite symmetric frameworks has also experienced some break- 
out results, building on a decade or more of initial Maxwell type necessary 



conditions for frameworks of various symmetry groups (Fowler and Guest 
2000} Guest and Fowler 2007} Connelly et al.}|2509) . In some key cases, these 



symmetry conditions predict finite motions for frameworks realized generi 
cally within the symmetry constraints, but whose graphs would be generically 
rigid without symmetry (Kangwai and Guest} [1999} |Bricard}p97) . Recently, 
key results of this work have been expressed in terms of 'orbit rigidity matrices' 
for the equivalence classes of vertices and edges under the group of symmetry 
operations S < |Schulze} [2009} |2010d{ |Schulze and Whiteley}[20Tu] . With modi- 
fied counts for the symmetry-preserving trivial motions ts, and with eo and vq 
denoting the number of edge orbits and vertex orbits under the group action 
of S, respectively, these matrices lead to Maxwell type necessary counts for 
frameworks to be infinitesimally rigid: eo ^ dvo - tg- 

Given that many crystal structures combine both periodic structure and 
symmetry within the unit cells, it is natural to investigate the interactions of 
these two types of group operations. So we will consider frameworks with 
'combined symmetry groups' of the form Z d x S, where Z rf is the group of 
translations of the framework, S is the group of additional symmetries of the 
framework, and x denotes the semi-direct product of S acting on Z d . Note that 
every symmetry operation in such a group can be written as a unique product of 
an element of Z d and an element of S. However, since S is typically not normal 
in Z d x S, the groups Z x S are in general not direct products. Details on the 
semi-direct product can be found in any abstract algebra text, such as Dummit 



and Foote ( |1991f . In Section[6]we will introduce combined 'orbit matrices' for 
the groups Z d x S. Combined with the counts of the trivial motions which 
preserve both symmetry and periodicity this will provide extended Maxwell 
type necessary counts for infinitesimal rigidity. In this setting we: 

1. count the rows of the combined orbit matrix: one row per orbit of edges 
r = e ; 

2. count the columns of the combined orbit matrix: one vector column per 
orbit of vertices plus columns for symmetry-preserving lattice deforma- 
tions: c = dvo + is', 

3. the dimension of the space of trivial motions (translations) left by sym- 
metries: is- 

The minimum dimension of the space of non-trivial symmetry-preserving in- 
finitesimal periodic motions of the periodic structure is: 

m = c — ts —r or m = dv® +€s~ts ~ e o- 

This is compared with the corresponding count on the graph without symmetry 
where with orbits of size k$ and no fixed edges or vertices, for the fully flexible 
lattice, we would anticipate: 

m = d(k s v Q ) + {^ + 2 J-d -(k s e ). 

In addition, if we choose the positions of the vertices generically within the 
symmetry (i.e., make one generic choice for each orbit of vertices) then the 



predicted infinitesimal motions will be finite flexes (Asimov and Roth 1978 
Schulze}[2uT0cT}|Schulze and WMeley}|2010) . 



The results are a surprise - adding symmetry can sometimes cause ad- 
ditional flexibility beyond what the original graph without symmetry would 
exhibit in the periodic lattice. These more flexible examples include symme- 
tries such as inversive symmetry, or half -turn symmetry with a mirror, found 
in a number of crystals, such as zeolites. Recent studies have confirmed that 
flexibility is a feature of natural zeolites (Kap kcTet al.| [2010 ) and contributes 
to their physical and chemical properties. In turn, this suggests that predicted 
flexibility in a computer designed theoretical 'zeolite' would be a criterion 
for selecting which theoretical compounds should be synthesized for further 
testing. 

When adding symmetry to a periodic lattice structure, we must consider 
the flexibility that this symmetry allows in the lattice structure. Inversive 
symmetry will be a key example, as it fits all possible lattice deformations (it 
occurs in 'triclinic lattices'), and the addition of this symmetry to the framework 
generates flexes from frameworks that previously were minimally rigid, while 
preserving the full range of possible flexes of the lattice itself. 
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In contrast, only certain types of lattices leave open the addition of a half- 
turn symmetry in 3-space. A half-turn parallel to a side of the lattice requires 
that side to be perpendicular to the remaining parallelogram face. This leaves 
only four of the six possible flexes of the lattice (monoclinic lattices), but it does 
predict additional flexes. Similarly, mirrors of symmetry can fit parallel to faces 
of the lattice, and restrict the shapes to monoclinic lattices, with the variable 
angle now parallel to the mirror. 

We can also have a larger symmetry group, with several generators. For 
example, monoclinic prismatic crystals, such as some forms of zeolite, have 
the symmetry group C2/1 which has both a half-turn symmetry and a mirror. 
These restrict the possible lattice shapes to lattices with a parallelogram base 
(perpendicular to the axis, parallel to the mirror) and the vertical prism at 
right angles to the base. Other forms of zeolite have the added symmetry of 
D 2 h, forcing the base to be a rectangle. Each symmetry group for the crystal 
structure and the associated crystal system requires some specific terms in the 
analysis, although patterns emerge, and we will present tables with rows for 
the combinations. 

In the larger theory of rigidity of frameworks, infinitesimal flexes of 'generic 
frameworks' transfer to finite flexes, for appropriate versions of generic. This 
holds for generic frameworks without symmetry, for frameworks generic 
within the symmetry class, and for periodic frameworks. That property ex- 
tends to these combined symmetry periodic frameworks, so we are talking 
about flexibility on a finite scale, at generic realizations for representatives of 
the orbits of the expanded group. 

The outline of the paper is as follows. In §2, §3, and §4, we present the basic 
definitions and associated rigidity matrices for: (i) finite frameworks (§2); (ii) 
finite frameworks with symmetry and the orbit matrices (§3); and (iii) periodic 
frameworks with the associated lattice matrices (§4). 

In §5, we introduce our method of inserting symmetry into the analysis of a 
periodic framework in dimensions 2 and 3. This includes a short summary of 
the wallpaper pattern types (dimension 2) and crystal systems (dimension 3) 
which arise as the possible lattices and restricted lattice variations for various 
symmetry groups in the corresponding dimension. §6 then presents some key 
examples in the plane for the groups Z 2 x C2 and Z 2 x C s with various lattice 
flexibilities, as well as summary tables over all plane groups of the form Z 2 x S. 
§7 presents key examples in 3-space, with corresponding tables. 

In §8 we briefly describe a range of extensions which are accessible using 
these methods, including: extensions to include fixed vertices and edges; exten- 
sions to additional plane and space groups; extensions to higher dimensions; 
and the companion static analysis of the frameworks. 

The analysis is not complete. While we have covered groups of the form 
Z d x S, some plane and space groups are not covered, namely those with 6- 
fold symmetry, or with glide reflections. As §8 illustrates, there is lots of room 
for additional exploration. We hope that this introduction, and the follow- 
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up papers with detailed proofs for what is claimed here, will offer interested 
researchers tools to explore the range of examples of interest in their context. 



2 Preliminaries on the rigidity of (finite) frame- 
works 

A framework in R d is a pair (G, p), where G is a finite simple (no loops or multiple 
edges) graph with vertex set V(G) and edge set E(G), and p : V(G) — > ]R d isamap 
such that p(z') + p(j) for all {i, j] e E(G). We also say that (G, p) is a d-dimensional 
realization of the underlying graph G (Whiteley J1996) . For i e V(G), we say that 



p(i) is the joint of (G,p) corresponding to i, and for {i,j} e E(G), we say 
that piy) := p({i,/}) is the bar of (G,p) corresponding to {i,j}. Moreover, we let 
v := |V(G)| and e := |£(G)|. It is often useful to identify p with a vector in 
by using the order on V(G). In this case we also refer to p as a configuration of v 
points in Mr. 

A framework (G,p) in R d is flexible if there exists a continuous path, called a 
finite flex or mechanism, p(f) : [0, 1] — > R' fo such that 

(i) p(0) = p; 

(ii) ||p(£); - p(t)j\\ = ||p, - Pj \\ for all < t < 1 and all e E(G); 

(iii) \\p(t)k - p(t)i\\ + \\pk - Pi\\ for all < t < 1 and some pair {fc,/} of vertices of 
G. 

Otherwise (G, p) is said to be rfgfd. For some alternate equivalent definitions of 
a rigid and flexible framework see Asimov and Roth ( |1978) , for example. 



An infinitesimal motion of a framework (G, p) in ]R rf is a function u : V(G) — > ]R d 
such that 

(p, - pj) ■ (m; - U; ) = for all ;} e £(G), (1) 

where m, denotes the vector w(f) for each i. 

An infinitesimal motion u of (G,p) is an infinitesimal rigid motion (or trivial 
infinitesimal motion) if there exists a skew-symmetric matrix S (a rotation) and a 
vector t (a translation) such that w, = Sp; + £ for all i e V(G). Otherwise u is an 
infinitesimal flex (or non-trivial infinitesimal motion) of (G,p). 

(G,p) is infinitesimally rigid if every infinitesimal motion of (G,p) is an in- 
finitesimal rigid motion. Otherwise (G,p) is said to be infinitesimally flexible 
( |Whiteley}[T9l6| |. 



While an infinitesimally rigid framework is always rigid, the converse does 
not hold in general. |Asimov and Roth| ( |1978[ l however, showed that for 'generic' 
configurations, infinitesimal rigidity and rigidity are in fact equivalent. 

The rigidity matrix of (G, p) (which in structural engineering is also known 



as the compatibility matrix of (G, p) ( jKangwai and Guest} |2000 Connelly et al 
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2009 1) is the e x dv matrix 



r(g, P ) ={;,/} 



(pi-pj) 



,j-pi) 



that is, for each edge \i, j) £ E(G), R(G, p) has the row with (p, -pj)\,..., (p, - p y )d 
in the columns d(i - 1) + 1, . . .,di, (pj - p;)i/ . . . , (pj - pi)d in the columns d(j - 1) + 
1, . . . , dj, and elsewhere (Whiteley 1996[ . See also Example 3.1. 

Note that if we identify an infinitesimal motion u of (G, p) with a column 
vector in R d " (by using the order on V(G)), then the equations in jl| can be 
written as R(G,p)w = 0. So, the kernel of the rigidity matrix R(G,p) is the space 
of all infinitesimal motions of (G, p). It is well known that the infinitesimal rigid 
motions arising from d translations and ( d ) rotations of R d form a basis of the 
space of infinitesimal rigid motions of (G,p), provided t hat the points m , ... ,p v 
span an affine subspace of R d of dimension at least d - 1 ( fWhiteley 1996} . Thus, 
for such a framework (G, p), we have nullity (R(G, pfj > d + (^) = (f^) and (G, p) 
is infinitesimally rigid if and only if nullity (R(G,p)) = C^ 1 ) or equivalently, 
rank(R(G,p)) = dv-{ d f). 

In particular, it follows that we can sometimes detect infinitesimal flexes in 
frameworks - and, by the result of Asimov and Roth ( 1978 ), even predict/i'nffe 
flexes in generic frameworks - by simply counting vertices and edges: 

Theorem 2.1 (Maxwell's rule) Let (G, p) be a d-dimensional framework whose joints 
span an ajfine subspace of R d of dimension at least d — 1. If 



e < dv - 



d + 1 



(2) 



then (G,p) has an infinitesimal flex. 

If the joints of(G,p) are in generic position, then there even exists a finite flex of 
(G,p). 



3 Symmetry in frameworks 

3.1 Symmetric frameworks and motions 

Given a finite simple graph G with vertex set V(G) = {!,... ,n), and a map 
p : V(G) — > ]R d , a symmetry operation of the framework (G, p) in R rf is an isometry 
s of R d such that for some a s e Aut(G), we have 

s(pd = Pa s (i) for all i e V(G), 



6 



where Aut(G) denotes the automorphism group of the graph G ([Hall} |1969 
Schulze 2010a bj. The set of all symmetry operations of a framework (G,p) 



forms a group under composition, called the point group of (G,p) (Bishop, 1973: 



Hall} 1969 Schulze 2010ajb i. Since translating a framework does not change its 
rigidity properties, we may assume wlog that the point group of a framework 
is always a symmetry group, i.e., a subgroup of the orthogonal group 0(R d ) 
( jSchulzellMualB} . 

Throughout this paper, we will highlight the Schoenflies notation for the 
symmetry operations and symmetry groups, as this is one of the standard 
notations in the l iterature (see |Bishop| (|1973); Connelly et al. (|2009); |Fowler 



|and Guest| j2000);|Guest and Fowler|(|2007) ; |Hall| ( |1969) ; |Kangwai and Guest 
( |1999| 12000) ; |Schulze| ( |2010a} |2009| |2010b) , for example). In the later tables for 



crystallographic groups, we will show three notations in parallel, to ensure 
clearer communication with multiple audiences. 

In the Schoenflies notation, the groups we will focus on in our examples and 
tables are denoted by C s , C„, Cnv, C„h, Ci, !D„, and T> nn - For dimension 2 and 3, 
C s is a symmetry group consisting of the identity Id and a single reflection s, and 
C n is a cyclic group generated by an n-fold rotation C„. The only other possible 
type of symmetry group in dimension 2 is the group C nv which is a dihedral 
group generated by a pair {C,„s}. In dimension 3, C nv denotes any symmetry 
group that is generated by a rotation C„ and a reflection s whose corresponding 
mirror contains the rotational axis of C„, whereas a symmetry group C„h is 
generated by a rotation C n and the reflection s whose corresponding mirror is 
perpendicular to the C„-axis. The group C, consists of the identity Id and an 
inversion i in 3-space. Finally, T) n is generated by an n-fold rotation C„ and a 
2-fold rotation C2 whose rotational axes are perpendicular to each other, and 
!D„h is generated by the generators C„ and C2 of a group D n and by a reflection 
s whose mirror is perpendicular to the C„-axis. 

Given a symmetry group S in dimension d and a graph G, we let &{Gg\ 
denote the set of all d-dimensional realizations of G whose point group is either 
equal to S or contains S as a subgroup (Schulze 2010a|b) . In other words, 
the set M(G£) consists of all frameworks (G,p) for which there exists a map 
cp:S^> Aut(G) so that 

s {pi) = P®(s)(i) f° r a ^ ! e all s e S. (3) 

If a framework (G, p) e &(g,s) satisfies the equations in |3| for the map O : S 



Aut(G), we say that (G,p) is of type O. It is shown in Schulze| ( [2009] |20 10b) that 
if the map p of a framework (G,p) € ^(G,S) is injective, then (G,p) is of a unique 
type O and <5 is necessarily also a homomorphism. For simplicity, we therefore 
assume that the map p of any framework (G,p) considered in this paper is 
injective (i.e., p, + pj if i + f). In particular, this allows us (with a slight abuse of 
notation) to use the terms p s in and po(s)(i) interchangeably, where i e V(G) and 
s e S. In general, if the type O is clear from the context, we often simply write 
s(z) instead of 0(s)(f). 
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An infinitesimal motion u of a framework (G, p) 6 &(g,S) is S-symmetric if 



s(m,) = « s (,) for all z e V(G) and all se5, 



(4) 



i.e., if u is unchanged under all symmetry operations in S (see also Figure [TJa) 
and(b)). 







(a) 



(b) 



(c) 



Figure 1: Infinitesimal motions of frameworks in the plane: (a) a C s -symmetric 
infinitesimal flex; (b) a C s -symmetric infinitesimal rigid motion; (c) an infinites- 
imal flex which is not C s -symmetric. 

Note that if for (G,p) e 3%<£S)i we choose a set of representatives {1, ...,Vo] 
for the orbits S(i) = {s(i)\s e S] of vertices of G under the group action of S, 
then the positions of all joints of (G, p) are uniquely determined by the positions 
of the joints pi, . . . ,p Va and the symmetry constraints imposed by S. Similarly, 
an »S-symmetric infinitesimal motion u of (G, p) is uniquely determined by the 
velocity vectors u\,..., u Vo for the representative vertices. 

The following extension of the theorem of |Asimov and Roth (1978 1 shows 
that an analysis of the '5-symmetric' infinitesimal rigidity properties of a sym- 
metric framework can be used to also detect finite flexes in the framework, 
provided that its joints are positioned generically within the symmetry. 



Theorem 3.1 ( Schulze 2010d > (see also Schulze ( |2009f ) Let S be a symmetry group 
in dimension d, and let (G,p) e ^(g,5) be a framework whose joints span all ofM. d , in 
an affine sense. If(G,p) is generic modulo the symmetry group S, i.e., the vertices of a 
set of representatives for the vertex orbits under the action ofS are placed in 'generic' 
positions (see Schulze §2010d\ty ; \Schulze and Whiteley ( 2010\ for details), and (G,p) 
also possesses an S-symmetric infinitesimal flex, then (G, p) also has a finite flex which 
preserves all the symmetries in S throughout the path. 



3.2 Orbit rigidity matrices for symmetric frameworks 



To determine whether a given framework (G,p) e &(G£) possesses an S- 
symmetric infinitesimal flex, we can use the techniques from group repre- 
sentation theory described in Fowler and Guest (2000|l; Kangwai and Guest 
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( 2000 ) ; [Schulze| ( |2010d|a) . In the recent paper Schulze and Whiteley| ( |20T0) the 
'orbit matrix' was introduced as a simplifying alternative to detect symmetric 
infinitesimal motions in symmetric frameworks and to predict finite flexes for 
configurations which are generic within the symmetry. In fact, it is shown in 



Schulze and Whiteley ( 2010 ) that t he orbit matrix is equivalent to the submatrix 



Schulze 



( 2010d I, but the construction is transparent, 
explicitly derived, without using techniques 



block Ri(G,p) studied in 
and the entries in the matrix can be t 
from group representation theory. 

For a d-dimensional framework (G,p) e &(GS) which has no joint that is 
'fixed' by a non-trivial symmetry operation in S (i.e., (G, p) has no joint p, with 
s(p,) = p, for some s e S, s + id), the construction of the orbit matrix becomes 
particularly easy (see Definition |3.1) , because, in this case, the orbit matrix has 
a set of d columns for each orbit of vertices under the group action of S. 



Definition 3.1 ( Schulze and Whiteley 2010 1 Let Shea symmetry group in dimension 
d and let (G, p) e £%(g$) be a framework which has no joint that is 'fixed' by a non-trivial 
symmetry operation in S. Further, let ffy(G) = {1, Col be a set of representatives for 
the orbits S{i) = {s(i)\s e S] of vertices ofG. For each edge orbit S(e) = {s(e)\ s e S] 
of G, the orbit matrix 0(G, p, S) of (G, p) has the following corresponding (dvQ- 
dimensional) row vector: 

Case 1: If the two end-vertices of the edge e lie in distinct vertex orbits, then there exists 
an edge in S(e) that is of the form {a,s(b)} for some s e S, where a,b e &v(G)- 
The row we write in 0(G, p, S) is: 



(0...0 (p a -s(p b )) 0...0 (p b -s-\p a )) 0...0). 

Case 2: If the two end-vertices of the edge e lie in the same vertex orbit, then there 
exists an edge in S{e) that is of the form {a, s(a)}for some s e S, where a e &v(G)- 
The row we write in 0(G, p, S) is: 

a 

(0...0 (2p fl -s(p fl )-s- 1 (p a )) 0...0). 



Example 3.2.1 To illustrate the above definition, we consider the 2-dimensional 
framework (G,p) with point group Ci — {id, C2) depicted in Figure [2] as an 
example. If we denote p\ = (a,b), p 2 = (c,d), p 3 = (-a,-b), and pi = (-c,-d), 
then the rigidity matrix of (G, p) is 



U,2} 
U,C 2 (2)} 
C 2 ({1,2}) 
C 2 ({1,C 2 (2)}) 



((a-c,b-d) 
(a + c,b + d) 
00 
00 



2 

(c - a,d - b) 
00 
00 

(a + c,b + d) 



3 = C 2 (l) 
00 
00 

(c — a,d — b) 
(-a — c,—b — d) 



4 = C 2 (2) 
00 

(-a -c,—b- d) 
(a - c,b - d) 
00 
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The orbit matrix 0(G,p,C2) of (G,p) will only have two rows, one for each 
representative of the edge orbits under the action of C 2 . (Note that if we are 
only interested in infinitesimal motions and self -stresses of (G, p) that are C2- 
symmetric, then it indeed suffices to focus on the first two rows of the rigidity 
matrix of (G,p). The other two rows are clearly redundant in this symmetric 
context!). Further, 0{G,p,C2) will have only four columns, because G has only 
two vertex orbits under the action of C2, represented by the vertices 1 and 2, 
for example, and each of the joints p\ and p2 has two degrees of freedom in the 
plane. Since both edge orbits satisfy Case 2 in Definition 3.1 0(G, p, C2) has the 
following form: 



1 2 
{1,2} ( (pi-p 2 ) (Pz-pi) ' 

{i,c 2 (2)} \(pi-c 2 {p 2 )) (pz-q 1 ^)) 



1 2 

(a — c,b — d) (c — a,d — b) 
(a + c,b + d) (c + a,d + b) 




(a) 




(b) 



Figure 2: The framework (G,p) € £%(g,c 2 ) ( a ) an d its corresponding symmetric 
orbit graph (b). 



We can use the 'symmetric orbit graph' to describe the underlying combina- 
torial structure for the orbit matrix of a symmetric framework (see also Figure 
0(b)): 

Definition 3.2 The symmetric orbit graph Gs of a framework (G,p) e &(g,S) is a 
labeled multigraph (it may contain loops and multiple edges) whose vertex set {l,...,v } 
is a set of representatives of the vertex orbits ofG under the action ofS, and whose edge 
set is defined as follows. For each edge orbit ofG under the action ofS, there exists one 



edge in G$: for an edge orbit satisfying Case 1 of Definition 3.1 Gs has a directed edge 
connecting the vertices a and b. If the edge is directed from a to b, it is labeled with s, 
and if the edge is directed from b to a, it is labeled with s~ l . For simplicity we omit the 
label and the direction of the edge ifs = id. Similarly, for an edge orbit satisfying Case 
2 of Definition \3.1\ Gs has a loop at the vertex a which is labeled with s. 

The key result for the orbit matrix is the following: 



Theorem 3.2 \Schulze and W hiteley, 2010) Let Sbea symmetry group and let (G,p) 
be a framework in &(G,sy Then the solutions to 0(G,p,S)u = are isomorphic to the 
space of S-symmetric infinitesimal motions of the original framework (G,p). 
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As an immediate consequence of Theorems 3.1 and 3.2 we have the fol- 
lowing Maxwell type counting rule for detecting finite 'symmetry-preserving' 
flexes in symmetric frameworks: 



Theorem 3.3 §Schulze\ fi009\ \2010d\ \Schulze and Whiteley\ |2010| Let S be a sym- 
metry group in dimension d and let (G, p) be a framework in S%(g$) which has no joint 
that is 'fixed' by a non-trivial symmetry operation in S. Further, let eo and Vo denote 
the number of edge orbits and vertex orbits under the action ofS, respectively, and let 
trivs denote the dimension of the space of S-symmetric infinitesimal rigid motions of 
(C,p). If 

e < dv - triv s , (5) 

then (G,p) has an S-symmetric infinitesimal flex. If the joints of(G,p) also span all of 
R d (in an affine sense) and are in generic position modulo S, then there even exists a 
finite flex of(G,p) which preserves the symmetries in S throughout the path. 

The dimension trivg of the space of »S-symmetric infinitesimal rigid motions 
of (G, p) e &( G,S) can easily be computed using the techniques described in 



Schulze ([2009 , 2010a). In particular, in dimension 2 and 3, trivs can be deduced 
immediately from the character tables given in |Connelly et~aL] p009) . Thus, 
in order to check condition |5}, it is only left to determine the size of the orbit 
matrix 0(G, p, S), which in turn requires only a simple count of the vertex orbits 
and edge orbits of the graph G under the action of S (see also Example 3.2.1 
and Table[l). 

Note that for a symmetry group S in any dimension d, the dimension ts of 
the space of S-symmetric infinitesimal translations can also be obtained in a very 
intuitive way, without using the techniques in Schulze (2009, 2010a I: if for a 



symmetry operation s e S, we let F s denote the symmetry element corresponding 
to s (i.e., F s = {a e R d \s(a) = a}), then t$ is simply the dimension of the symmetry 
element of the group S, i.e., the dimension of the linear subspace OseS Fs of JR d / 
because the initial velocity vectors of an 5-symmetric infinitesimal translation 
must all be contained in the space HlseS F s . 

For example, for the 'reflectional' symmetry group C s = {id, a} in dimension 
d, it is easy to see that the space of C s -symmetric infinitesimal translations is of 
dimension (d — 1), since it consists of those translations whose velocity vectors 
are elements of the (d - l)-dimensional mirror-plane F a - R d n F a = Fid H F a = 
Hsec, F s corresponding to a (see also Figure [T](b)). 

However, finding the dimension of the space of 5-symmetric infinitesimal 
rotations heuristically without the techniques in Schulze] < |2009| p010a I, becomes 
increasingly hard, if not impossible, in dimensions > 3. 



Example 3.2.2 Let's apply Theorem |3.3| to the framework (G, p) we considered in 
Example 3.2.1 (see also Figures|2](a) and|3j. We clearly have dvo =2-2=4 and 
eo = 2. Further, we have trivc 2 = 1, since the only infinitesimal rigid motions 
that are C2-symmetric are the ones that correspond to rotations about the origin 
(see Schulze < |2009 2010a} for details). Thus, we have 

e = 2 < 3 = dv - trivc 2 ■ 
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(a) (b) 



Figure 3: A C2-symmetric infinitesimal flex of the framework from Example 
3.2.1 (a) and the path taken by the joints of the framework under the corre- 
sponding symmetry-preserving finite flex (b). 



So, by Theorem |3.3| we may conclude that any realization of G which is 'generic' 
modulo the half-turn symmetry has a symmetry-preserving finite flex (Fig- 
ure^. □ 
Note that the standard (non-symmetric) Maxwell count also detects a finite 
flex for the framework in Example 3.2, since for the graph G, we have e = 4 < 5 = 
2v - 3. However, as shown in Schulze and Whiteley (2010), there exist a range 



of interesting and famous examples in 3-space, including the Bricard octahedra 
< |Bricard[ |Stachel} [1987) ) and the flexible cross-polytopes ( |Stachel[ [2000) , 
which can be shown to be flexible via the symmetric count in Theorem |3.3| but 
not with the standard non-symmetric Maxwell-Laman type counts for rigidity. 

The following table shows the symmetric Maxwell type counts for a selec- 
tion of point groups S in 3-space. For simplicity at this stage, we assume that 
no joint and no bar is fixed by a non-trivial element in S, so that all vertex orbits 
and edge orbits under the action of S have the same size kg. (Recall that a joint 
Pi is fixed by s e S if s(p,) = p,; a bar {p;,p ; } is fixed by s e S if either s(p,) = p, 
and s(p ; ) = pj or s(p,) = pj and s{p j) = pi). So, in particular, both the number of 
joints, v, and the number of bars, e are divisible by kg. A necessary condition 



for rigidity in 3-space is e > 3c - 6 (recall Theorem 2.1 1. For each group S 
in Table [l] e is chosen to be the smallest number which satisfies e > 3v - 6 
and is divisible by k$; that is, e is chosen to be the least number of edges for 
the framework to be rigid without symmetry and to be compatible with the 
symmetry constraints given by S. The integer fg in the final column indicates 
an fg -dimensional space of 5-symmetric infinitesimal flexes if fg > 0, and a 



(-/s)-dimensional space of »S-symmetric self-stresses if fs < 0. (See Section 8.5 
for more on self-stresses.) 

The final column of Table [T] indicates that at 'generic' configurations, the 
frameworks with Ci symmetry always have a finite flex, while those with C4 
symmetry are always stressed. 



12 



Table 1: Impact of some 3-space point groups on counts for rigidity. 



s 


ks 


trivs 


e 


e o 


3v - 


trivs 


fs 


Ci 


1 


6 


3v-6 


3v - 6 


3v 


-6 





Ct 


2 


3 


3v-6 


3v - 3 


3v 


-3 





c 2 


2 


2 


3v-6 


3v - 3 


3v 


-2 


1 


C s 


2 


3 


3v-6 


3v - 3 


3v 


-3 





C-2M 


4 


1 


3v-A 


3v - 1 


3v 


- 1 





®in 


8 





3v 


3v 


3v 


-0 







4 


2 


3v-A 


3v - 1 


3v 


-2 


(-1) 


c 3 


3 


2 


3v-6 


3^o-2 


3v 


-2 






4 Periodic frameworks 

4.1 From infinite periodic frameworks to a finite orbit graph 

Let G be a simple infinite graph with finite degree at each vertex. Let p be a 
placement of the vertices V(G) in ]R 3 , such that the resulting framework (G, p) is 
invariant with respect to three linearly independent translations ti,t2,t3 e R 3 . 
We assume without loss of generality that ti lies on the x-axis, and t2 lies in the 
xy-plane. Let L be the matrix whose rows are these translations: 





ti 




" t n 








L = 


t 2 






hi 







t 3 




hi 


hi 


f.33 



We call the pair ((G,L),p) a periodic framework. This definition can be adapted 
to describe periodic objects in d dimensions, but here our focus will be on 
crystal-like structures in two and three dimensions. In two dimensions, 



L = 



h 




" hi 





t 2 




hi 


h 2 



and the other definitions are similarly adapted. 

The three translations ti, t2, t 3 define a parallelepiped called the unit cell. This 
can be equivalently described by three lengths a, b, c, and three angles a, f>, y, 
as illustrated in Figure [4] These coordinates are standard in crystallography 
texts (Wikipedia| |2010[ |. We may obtain one representation from the other by a 
change of coordinates. 

A copy of the unit cell centered at the origin is given by 

U = (fliti + « 2 t2 + a 3 t 3 | - \ < at < i) , 
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Figure 4: The unit cell defined by the translations ti,t2,t3. The unit cube (all 
side lengths 1) is shown in the dashed lines. 



with its boundary defined to be the boundary of the closed parallelepiped: 

U = jfliti + a 2 t 2 + a 3 t 3 | |flf| < ^ j . 

We assume without loss of generality that no vertex of V(G) lies on the boundary 
of the unit cell (we may simply translate the framework ((G, L), p) until no vertex 
lies on the boundary). 

The translations ti,t2,t3 generate a crystalline lattice t]Z X t2Z X t 3 Z. This 
lattice partitions R 3 into copies of the unit cell, each containing exactly one 
lattice point. Let each cell be centered at that lattice point, and let the cells be 
indexed by Z 3 according to the lattice point they contain. 

We now use this partition to define a finite labeled graph G which represents 
the periodic framework, and can be used to study its rigidity. The edges of this 
graph are labeled invertibly by elements of thegroup Z 3 in a way that captures 
the periodic structure of ((G,L),p). Figure pi illustrates this process for the 
analogous 2-dimensional case, in which edges are labeled by elements of the 
group Z 2 . 

Definition 4.1 The periodic orbit graph (G,n) of a periodic framework {{G,L),p) 
is a labeled multigraph whose vertex set V(G) = {!,..., v] consists of the vertices 
appearing in the (0,0,0)-cell of ((G,L),p). The edge set E(G) and the labelling n : 
E(G) — » Z 3 on these edges are defined as follows: 

Case 1: An edge~ee E(G) whose length is completely contained within the (0, 0, 0)-cell 
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(a) «G,L>,p) (b) (G,n) 



Figure 5: A 2-dimensional periodic framework (a) with the (0, 0)-cell indicated. 
The periodic orbit graph is shown in (b). 



must connect two distinct vertices ofV(G) which also lie in the (0, 0, 0)-cell. For 
every such edge ~e, define e e E(G) to be the edge connecting the corresponding 
vertices ofV(G). Assign this edge an arbitrary direction, and label it with (0, 0, 0). 
For visual simplicity in our diagrams, edges labeled by the zero (identity) element 
appear as unlabeled, undirected edges. 

Case 2: Let ~e e E(G) be an edge that crosses the boundary of the (0, 0, 0)-cell. In 
particular, suppose T connects the vertex v in the (0,0,0)-cell with the vertex 
w in the (n\, ni, n?)-cell. Then define e e £(G) to be a directed edge of G that 
originates in the vertex v, terminates at the vertex w. Assign this directed edge 
the label (n\,ni,n^ 6 Z 3 . 



The periodic orbit graph (G,n) contains a set of representatives of the vertex 
and edge orbits of ((G,L),p) under the action of Z 3 . Let p be the restriction of 
p^to the vertices of the (0, 0, 0)-cell, p : V -» U. We call the pair ({G,n),p) the 
periodic orbit framework of ((G,L),p), and the labeled multigraph (G,n) will be 
called the periodic orbit graph. 

Let the edges E(G) be ordered. An edge of (G,n) is denoted by {i,j;n^, 
where i, j e V(G) and e Z 3 . This edge corresponds to an equivalence class of 
bars in the periodic framework ((G, L),p), which contains the bar (pi,Pj + n^T), 
and all of its translates by a,t„ for integers This edge can be equivalently 
represented by {/, i; -n^}. 

It has been shown (Malestein and Theran 2010 Ross 20101 that every (d- 
dimensional) periodic graph admits such a representation, and is invariant 
under the choice of unit cell of a specified size. It is also known that every 
finite directed multigraph whose edges are labeled wit h elements of the group 
Z d can be realized as a rf-dimensional periodic graph I Ross 2010k. In general, 
directed multigraphs whose edges are labeled by elements of a group, with the 
reverse direction implicitly labeled by the inverse group element, are known as 
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gain graphs (Gross and Tucker 



2001) . 



4.2 Periodic rigidity and infinitesimal rigidity 

We may define notions of rigidity and infinitesimal rigidity for the periodic 
framework ((G, L),p) simply by a direct application of the ideas from Section]^ 
However, in this paper we are interested exclusively in motions and infinitesi- 
mal motions of ((G, L), p) which preserve the periodicity of the framework. We 
say that ((G,L),p) is periodic (infinitesimally) rigid if the only periodic (infinites- 
imal) motions of ((G, L),p) are trivial. By defining rigidity and infinitesimal 
rigidity for a periodic orbit framework ((G,n),p), we are able to identify pre- 
cisely these characteristics of the corresponding periodic framework, «G, L), p). 
Further details can be found in |Borcea and Streinu] ( |2010bf ; Malestein and| 
Theran| ( |20T0) ; [Rcis1 ( |20T0) . 



Let ((G, n), p) be a periodic orbit framework, and let e k = [i, j; n k } be an edge 
of (G, n), with n k = (n; cl , n k2 , %). The edge length of e k is given by the Euclidean 
length of the vector p,- - (pj + n k L). 

3 

Ifoll 2 = lip, - (Pi + n k L)\\ 2 = Y^iPre - (Pjc + (n k L) f )) 2 

e=\ 

= (Pa - (P;i + n kl hi + nnh\ + n^hi)) 2 + 

(pc - (P/2 + n k2 t 2 i + n k3 t 32 )) 2 + (p, 3 - (p/3 + n k3 t 33 )) 2 (6) 

Because there are a finite number of edges of (G,n), we may use the above 
periodic edge length to define rigidity and flexibility of the orbit framework 
((G,n),p). It is analogous to the definition presented in Section|2]and we omit 
it here. 

Letting the positions of the vertices p, = (pn,pi2,pi3) and the generators of 
the lattice t\\, . . . , t 33 vary with time, infinitesimal motions of ((G, n), p) can be 
found by differentiating j6|, with the assumption that ||e/ c || = C, a constant. We 
obtain 

(p, - (pj + n k L)) ■ (dpi/dt - dpj/dt) + £(i, j; n k ) ■ (dt n /dt, dt 33 /dt) = 0, 

where £,(i, j; n k ) is the 6-tuple of coefficients of (dt\\ /dt, ... , dt 33 /dt) . For example, 
the coefficient corresponding to t 32 is n k3 (pa - (pj 2 + n k2 t 22 + n k3 t 32 )). 

More generally, an infinitesimal motion (u, w) of a periodic orbit framework 
«G, n), p) in R 3 is a pair of functions 

u : V{G) M 3 , and w : L -> R 6 

such that 

(pi - (pj + n k L)) ■ (ui - Uj) + £(i, j; n k ) ■ w = for all [i, j; n k \ 6 E((G, n». (7) 
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An infinitesimal motion (u,w) of ((G,n),p) is called trivial if w is the zero map, 
and there exists uq 6 M 3 such that u(i) = uq for all i 6 V(G). This simple form 
follows from the fact that the only isometries of the whole space that preserve 
the periodic structure are translations. If an infinitesimal motion is not trivial, 
then it is called an infinitesimal flex. The periodic orbit framework ((G,n),p) 
is infinitesimally rigid if every infinitesimal motion of ((G, n), p) is a trivial one. 
Otherwise the framework is infinitesimally flexible. 

Let (u,w) be an infinitesimal motion of ((G,n),p), with u + 0. When w : 
L — > R 6 is the zero-map, it indicates that the lattice vectors are fixed by the 
infinitesimal motion. We call such a motion («,0) a lattice-fixing infinitesimal 
motion. If w + 0, then the infinitesimal motion (u,w) is called lattice flexing. In 
this paper we regard the lattice-fixing motions as a specialization of the lattice- 
flexing motions, and hence we will assume all motions are lattice-flexing, unless 
otherwise noted. 

The fixed lattice variation is interesting in its own right, and in two- 
dimensions admits a concise combinatorial characterization (Ross |20Tu) . It 



may also be of interest to consider partial flexing of the lattice. For example, 
we may ask that the three translation vectors are only allowed to scale, but 
the angles between them remain fixed at 90°. This would correspond to the 
translation matrix 

£n 
£ 2 2 
£ 33 

Such variations will be treated briefly at the end of this section. 

Remark 4.1 Any infinitesimal motion of a periodic orbit framework ((G,n),p) can 
be extended to a periodic infinitesimal motion of the periodic framework ((G,L),p). 
However, there will be some infinitesimal motions of((G, L),p) that do not preserve the 
periodic structure, and therefore do not specialize to infinitesimal motions of({G, n), p). 
In particular, any infinitesimal motion of ((G,L),p) that breaks the periodicity of 
((G,L),p) will not appear as a motion of ((G,n),p). In other words, a periodic frame- 
work ({G, L),p) may be infinitesimally flexible, yet infinitesimally periodic rigid in our 
analysis. 



It has been demonstrated Borcea and Streinu (2010b I; Malestein and Theran 



< |2010[ >; Ross < |2010 ) that the vector space of periodic infinitesimal motions of a 
periodic framework {{G, L),p) is equivalent to the vector space of infinitesimal 
motions of the periodic orbit framework {{G,n),p). To make this connection, 
we use a periodic orbit matrix, which we shall now define. 

4.3 Orbit rigidity matrices for periodic frameworks 

Let {{G,n),p) be a three-dimensional periodic orbit framework with v = \V{G)\ 
and e - \E(G)\. The rigidity matrix R((G, n),p) for the periodic orbit framework 



17 



is the e x (3c + 6) matrix with the row corresponding to the edge {i, j; n k \ given 
by 

i l 

U,j;n k } (0---0 (pi-(p } + n k L)) 0---0 ( P/ - (p { - n k L)) 0---0 £(i,j;n k )). 

The entry £,{i,j;n k ) is a six-tuple representing the coefficients of 
(dtn/dt, . . . ,dt^/dt). Loops may appear in the periodic orbit framework. The 
row of R((G, n), p) corresponding to the loop edge {i, i; ne) is 

i L 
(0---0 0---0 0---0 £(i,i;n { )). 

Note that the matrix R((G,n),p) is identical to the 'augmented compatibility 
matrix' used by |Guest and Hutchinson ( |2003[ |. 



Example 4.3.1 We consider a two-dimensional periodic orbit framework with 
two vertices and five edges. The rigidity matrix R((G, n),p) for this framework 




(0,2) 




(1,1) (1,0) 



(a) {(G,L),p) (b) (G,n) 

Figure 6: A two vertex example in R 2 . 



has 5 rows and 7 columns. Below it is broken into two sections: the four 
columns corresponding to variables p\ = (a, b), pi = (c, d), and the three columns 
corresponding to the three non-zero variables in L: fn,f2i,f22- Note that the 
edge {1, 2; (0, 0)} has the entry (0, 0, 0) in the (fu, tzi, hi) columns, but is non-zero 
elsewhere. In contrast, the loop edge {1,1; (1,0)} has zero entries everywhere 
except in the columns corresponding to L. 





Pi = (a,h) 


Pi = (c,d) 


(hi, hi, hi) 


{1,2; (0,0)} 


(pi " Pi) 


(pi ~ Pi) 


(0,0,0)) N 


{1,2; (0,2)} 


(pi-[p 2 + (0,2)L]) 


(p 2 -[pi-(0,2)L]) 


(*, *, *) 


R«G,n),p) ={2,1; (1,0)} 


(pi-[p 2 -(l,0)L]) 


(p 2 -[p 1 + (l,0)L]) 


(*, *, *) 


{2,1; (1,1)} 


(pi-[p 2 -(l,l)L]) 


(p 2 -[p 1 + (l,l)L]) 


(*, *, *) 


{1,1; (0,1)} 


I o 





(0,0,1) J 
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where the three columns corresponding to (tn, hi, £22) (the non-zero entries of 
L) are given by: 



hi 



-(c - (a + tn)) 

-11 + ^21)) 



-(c - (a + 1 




-(c- 



hi 




(a + t n + £21)) 




£22 


-2(b - (d + 2*22)) 


~(d ~(b + £22)) 
1 



The two trivial motions of ((G,n),p) are represented by the column vectors 
(1, 0, 1, 0, 0, 0, 0) r and (0, 1, 0, 1, 0, 0, 0) T , which are always solutions to the linear 
system R((G,n),p) • (u,w) = 0. These are translations of the whole structure. 
If we are only interested in the lattice-fixing infinitesimal motions of ((G, n), p), 
we may omit the columns corresponding to (tn, hi, hi) and the translations still 
appear in the modified matrix. □ 

Returning to three dimensions, we may associate an infinitesimal motion 
(u, w) of the periodic orbit framework ((G, n), p) with a column vector (u, w) T in 
]R 3d+6 . The equations in |7| may then be written as the solutions to the linear 
system described by R((G, n),p) • (u, w) T = 0. There will always be three trivial 
solutions corresponding to the three trivial motions, and hence the maximum 
rank of R((G,n),p) is 3c + 3. 

If we wish to consider only the lattice-fixing infinitesimal motions of 
((G,n),p), then we may omit the final six columns corresponding to £,(i,j;rik). 
Our rigidity matrix is then of dimension e X 3c, with maximum rank 3c - 3. In 
general, the maximal rank of a fixed-lattice rigidity matrix is dv — d, and the 
maximal rank of a flexible-lattice rigidity matrix is dv 
Streinu||2010bi. 



d + ( d f) (Borcea and 



As in the symmetric setting, there also exists a modified notion of generic 
for periodic frameworks (for further details see Ross (2010 1). For our purposes, 
it will be important to know only that generic rigidity of the periodic orbit 
framework ({G,n),p) depends on the underlying periodic orbit graph (G,n). 
Furthermore, for generic frameworks, infinitesimal rigidity and rigidity are 
equivalent (Malestein and Theran 2010 Ross |2010) . 

As shown in |Borcea"a nd Streinu (2010b I; Males tem and Theran| ( |2010[ |; |Ross 
( |2010| >, the vector space of periodic infinitesimal motions of a periodic frame- 
work ((G, L), p) is equivalent to the vector space of infinitesimal motions of the 
periodic orbit framework ((G, n), p). The following theorem states that the vec- 
tor space of periodic infinitesimal motions of a periodic framework corresponds 
to the kernel of the periodic rigidity matrix. 



Theorem 4.1 ^Borcea and Streinu 2010b Ross\ 2010\ Let ({G,n),p) be a periodic 
orbit framework corresponding to the periodic framework ((G,L),p). The kernel of 
the corresponding periodic rigidity matrix R({G,n),p) is isomorphic to the space of 
periodic infinitesimal motions of the associated periodic framework ((G,L),p). 
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Corollary 4.2 ( Borcea and Streinu 2010b) The periodic framework ((G,L),p) is in- 

4iMi1-0eiwi nll-t i v\ovir\Air \'1 m A 1 1 -j TT? ^ i-f nn A /i -n 7i / if ±Vto vni/i V- r\f 4-tio vim Ai4-1 1 Mi/rfi'i v £i v + 



finitesimally periodic rigid in W if and only if the rank of the rigidity matrix for the 
corresponding orbit framework R((G,n),p) is dv — d + C^ 1 )- 

Returning to Example 4.1, the rank of the matrix corresponding to generic 
positions of the vertices is 5, which is maximal on two vertices, and hence 
((G, n), p) is infinitesimally rigid. If we are only interested in the rigidity of the 
framework on a fixed lattice, the rank of the lattice-fixing portion of R((G, n),p) 
(the first four columns) is 2, which again is maximal. Note that this means that 
three of the five edges of our example are redundant on a fixed lattice. 



From Theorem 4.1 and a periodic version of Theorem|3.1|we obtain a peri- 



odic Maxwell type counting rule for detecting finite periodic flexes: 

Theorem 4.3 Let ((G,L),~p) be a periodic framework in dimension d with a corre- 
sponding orbit framework ((G,n),p), where v = \V(G)\ and e = \E(G)\. If 

e<dv — d + fTl = dv- 



2 ) \2)' 

then ({G,n),p) has an infinitesimal flex, which corresponds to a periodic infinitesimal 
flexof{{G,L),p). 

Furthermore, for generic positions of the vertices of G relative to the generating 
lattice L, ((G,n),p) has a finite flex, which corresponds to a periodic finite flex of 
«G,L),p). 



Theorem 4.3 can be adapted for the fixed lattice with the count e < dv - d, 
or for any other variation of the flexible lattice. For dimensions 2 and 3, Table 
[2] shows the number of lattice parameters corresponding to each of the lattice 
variants in the following list: 

(i) fully flexible lattice: all variations of the lattice shape are permitted; 

(ii) distortional change: keep the volume fixed but allow the shape of the lattice 
to change; 

(iii) scaling change: keep the angles fixed but allow the scale of the translations 
to change independently; 

(iv) hydrostatic change: keep the shape of the lattice unchanged but scale to 
change the volume; 

(v) fixed lattice: allow no change in the lattice. 

Why might we study one of these variants? For crystals, we might focus 
on short time-scale vibrations, during which the large motions of distant atoms 
needed for a flexible lattice could not happen. In this case we effectively study 
a fixed lattice with local variation and all velocities small. Or we might study 
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Table 2: Number of parameters corresponding to types of lattice deformations 
with no added symmetry, in two and three dimensions. 



LatticeDef 


2-D 


3-D 


flexible 


3 


6 


distortional 


2 


5 


scaling 


2 


3 


hydrostatic 


1 


1 


fixed 









slow responses to general pressure, given by a fully flexible lattice. In be- 
tween, we could study responses to pressures and constraints of various types, 
with various boundary conditions, which correspond to various intermediate 
situations. 

Another setting which produces periodic structures is simulation of large 
sphere packings by simulations with a modest number of spheres, and a peri- 
odic bounding box to give a better approximation than a fixed boundary. Here, 
which case applies will depend on the variation of the periodic bounding box 
which the simulation chooses to permit. 



5 Periodic frameworks with symmetry 

In the previous two sections we have built up the orbit matrix for finite frame- 
works under point groups S in R d and for periodic frameworks with groups 
Z d . Counting the rows and columns of these orbit matrices, the Maxwell type 
counts of rows vs columns minus trivial motions give necessary conditions for 
a framework to be generically rigid (minimally rigid). Recall that for Z d we 
had several variants, ranging from the fully flexible lattice with C^ 1 ) columns 
added for the lattice variables, to the fixed lattice with no columns added for 
lattice deformations. 

We now turn our attention to periodic frameworks with added symmetry. 
These frameworks have orbit graphs whose edges are labeled by elements 
of groups of the form Z d x S. An example of such a framework with its 
corresponding orbit graph is shown in Figure [7] 

Definition 5.1 Let ({G, L),p) be a periodic framework with symmetry group Z d x S. 
The symmetric periodic orbit graph (G,g) corresponding to this framework is the 
labelled multigraph with one representative for each equivalence class of edges and 
vertices under the action ofZ? x S. The labelling of the edges g : £(G) — > Z d x S is 



determined in the manner described in Definitions 3.2 andvLl 



The symmetry group Z d x S of a symmetric periodic framework will de- 
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Figure 7: A plane framework with Z 2 x Ci symmetry can be labeled with the 
elements of the group (a), or in short hand with 'gains' (b) as in the gain graph 
(c) 



termine its crystal system, which is a characterization of the parameters which 
determine the unit cell. This is usually defined by the number and arrangement 
of lengths and angles determining the cell, and these parameters represent the 
variations of lattice shapes which preserve the given symmetry ( |Wikipecfla 



2010} International Union of Crystallography , 2006). In the plane we will con- 



sider four different crystal systems, shown in Figure|8} and in space we consider 
the six crystal systems shown in Figure [9] 








b 
"1 






(a) square 



(b) rectangle 



(c) rhombus 



(d) parallelogram 



Figure 8: The four planar crystal systems. The number of lattice parameters 
are (a) 1, (b) 2, (c) 2, (d) 3. 
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fa) cubic 



(b) tetragonal 



(c) orthorhombic 




(d) trigonal (e) monoclinic (f) triclinic 



Figure 9: The six crystal systems addressed in this work. The number of lattice 
parameters are (a) 1, (b) 2, (c) 3, (d) 2, (e) 4, (f) 6. 

The crystal system of a framework specifies the maximum number of param- 
eters that determine the lattice, and therefore the number of lattice columns of 
our orbit rigidity matrix. We may further reduce the number of lattice columns 
by changing the type of lattice that we are considering: (flexible, distortional, 
scaling, hydrostatic or fixed), although it should be noted that the lattice sys- 
tem will partially determine these choices. For instance, for a two-dimensional 
framework with a rhombus unit cell, scaling and hydrostatic will be identical. 

Remark 5.1 If we transform the lattice to the unit cube, by an affine transformation 
which preserves the symmetry, then the lattice parameters represent the number of 
non-zero partial derivatives of the length of bars, for variables from the lattice pattern. 

The rigidity of the symmetric periodic framework can be studied using the 
orbit matrix. Letting Vq and eo represent the number of vertices and edges in the 
orbit graph G, the orbit matrix has dimension eo X (dvo + where ts describes 
the number of columns corresponding to the lattice parameters. This number 
will vary depending on a) the crystal system corresponding to the symmetry 
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group S, and b) the type of lattice variation we are considering (i.e. flexible, 
distortional, scaling, hydrostatic, or fixed). 

As with the orbit matrices for S, the second number which is important is 
the number of trivial motions which preserve all the group operations. Since 
we are working with periodic frameworks, we are looking for what translations 
also preserve the symmetries in S within the orbit matrix. We denote by t$ the 
dimension of the space of points which is fixed by all elements of the group 
S. This space is also called the symmetry element of the group. For our 
calculations here, this can only be one point, a line, a plane, or all of 3-space. 

It is the combination of these two numbers Is an d tsi plus the number 
of orbits of edges and vertices (corresponding to the number of edges eo and 
vertices Vq of the orbit graph (G,g)), which generates the predictions of the 
number of non-trivial motions (if any) which occur when the vertices of the 
framework are in generic position. As before, we assume that no edges or 
vertices are fixed by the action of the group Z rf x S. 

In Section[6]we describe symmetric periodic frameworks with two samples 
of S in the plane in detail, followed by tables covering all the groups within our 
analysis. In Section [7] we outline two samples of S in 3-space, again followed 
by tables for the relevant groups. 

In the tables we give three distinct notations for each group: the Schoenflies 
notation used by chemists ( jBishop 1973 Hall} 1969(1, the Hermann-Mauguinno 
tation used internationally by crystallographers ((International Union of Crys 



tallography 2006} , and the orbifold notation used by mathematicians (Conway 
et al-T|200ip . 



Recall that the Schoenflies notation was briefly introduced in Section 3.1 



In the Hermann-Mauguin notation, n-fold rotational symmetries are de- 
noted by n, and these axis numbers are written down in decreasing order of n. 
If an n-fold rotational axis is contained in a mirror, then an m is written after 
the corresponding axis number. If there exists a mirror which is perpendicular 
to an n-fold rotational axis, then this is denoted by placing the symbol jm after 
the corresponding axis number. Finally, the notation n for an axis number n 
indicates that an inversion, followed by the n-fold rotation, is a symmetry of 
the structure. 

In the orbifold notation for wallpaper groups in the plane, an n-fold ro- 
tational symmetry is denoted by n, and a mirror symmetry is denoted by an 
asterisk, *. A number before an asterisk indicates a center of pure rotation, 
whereas a number after an asterisk indicates a center of rotation with mirrors 
through it. If there are only translational symmetries present, then this is de- 
noted by the symbol o. The orbifold notation for groups in 3-space is similar; 
note, however, that for a group in 3-space, an x indicates the presence of an 
inversion symmetry, whereas in 2-space, x indicates a glide reflection. 

For an interdisciplinary audience, the simultaneous use of these three nota- 
tions seems to be an appropriate presentation, and such comparative columns 
can be found in multiple sources, including the Wikipedia pages for crystal 
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systems ( Wikipedia, 2010). 

Note that the assumption that the group has the form Z d x S means that 
we do not cover all the plane wallpaper groups, or the full set of space groups. 
Specifically, we will not include groups which have glide reflections as genera- 
tors of the group, or those which have 6-fold rotations. These other groups will 
also have orbit matrices, but require an alternative analysis for comparisons 
and counts. We return to this issue in Section|82] 



6 2-D periodic frameworks with symmetry: Z/ 2 x S 

6.1 Z 2 x C2 - half-turn symmetry in the plane lattice 

Half -turn symmetry in the plane is equivalent to inversion in the point axis. 
This symmetry fits an arbitrary parallelogram for the lattice (Figure |8|d)), and 
£c 2 = 3. We will consider periodic plane frameworks with symmetry Z 2 x C2 
for two variations of the lattice: (1) a fully flexible lattice; (2) a fixed lattice. 

Example 6.1.1: fully flexible lattice Z 2 x C2. The original (non-symmetric) 
necessary count for a periodic framework on the fully flexible lattice to be 
minimally rigid is e = 2v + l (recall Theorem fP) . To permit half -turn symmetry, 
with no vertex or edge fixed by the half-turn, we will need to start with the 
modified count 2eo = 2(2co) + 2, where vq and eo are the numbers of vertices and 
edges of the orbit graph, respectively. Dividing by 2, this gives eo = 2c>o + 1- 

Under the half -turn symmetry with a fully flexible lattice, the orbit matrix 
has 2 columns under each orbit of vertices, plus (c 2 = 3 columns for the three 
parameters for the lattice deformations. Further, we clearly have fc 2 = since 
there are no infinitesimal trivial motions which preserve the half -turn symmetry 
along with the periodic lattice. This creates the necessary symmetric Maxwell 
condition 

e > 2v + 3 - 

for periodic rigidity. However, as shown above, for a graph that was previously 
minimally rigid without the symmetry, we have eo = 2z>o + 1 < 2co +3. This gap 
predicts that a graph which counted to be minimally rigid without symmetry, 
realized generically with half-turn symmetry on a fully flexible lattice, now 
has two degrees of (finite) flexibility. As an example, consider the snapshots 
of three configurations with the same edge lengths but changing angles and 



lengths of the unit cell in Figure 10 Together these snap shots confirm the 
predicted two degrees of freedom. 



The orbit matrix corresponding to the framework pictured in Figure 10 has 
the following form: 
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Figure 10: A generically rigid graph on a fully flexible lattice, realized with 
2-fold symmetry has several non-trivial flexes changing the lattice. Its periodic 
symmetric orbit graph is pictured in (d). 





Pi 


P2 


P3 


(^11,^21,^22) 


{1,2} 


(pi - Pz) 


(P2 - Pl) 





(0,0,0) 


{2,3} 





(P2 - Pi) 


(P3 - P2) 


(0,0,0) 


{3,1} 


(pi " Ps) 





(P3 " Pl) 


(0,0,0) 


{1,2; gi} 


(pi-gityzj) 


[Pl-g^ipi)) 





(*, *, *) 


{1,2; g 2 ) 


(pi - giiptj) 


[pi-g^iPl)) 





(*, *, *) 


{2,l;g 3 } 


(pi - g^ipii) 


(P2 - g$(Pl)) 





(*/ *, *) 


{3,l;g 4 } 


Xpi-gl l (pz)) 





(P3 - gliPlj) 


(*, */ *) 



□ 

Example 6.1.2: fixed lattice Z 2 x Ci- The original (non-symmetric) necessary 
count for any minimally rigid periodic framework on the fixed lattice is e = 
2v - 2. With added C2 symmetry, we have e = 2eo and v = 2vo (all orbits have 
kc 2 = 2), so a minimally rigid orbit graph, realized with C2 symmetry, will have 
le = 2(2» )-2,ore =2» -l. 

As the example below illustrates, with a fixed lattice, the orbit matrix has 2 
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columns under each orbit of vertices. Further, there are no translations which 
preserve the half-turn symmetry along with the periodic lattice, and hence we 
have tc 2 = 0. This creates the necessary symmetric Maxwell condition 

e > 2v Q 

for periodic rigidity. However, as shown above, if the graph was chosen to be 
minimally rigid without the symmetry, we have eo = 2co - 1. The gap eg = 
2vo - 1 < 2co shows that with the added half -turn symmetry, a minimally rigid 
graph will become flexible within the fixed lattice. Figure [TT| shows the sample 
graph already presented in Figure [7] with vq = 3, eg - 5 and eo = 5 < 6 = 2vo, 
with two realizations with the same edge lengths - illustrating snapshots of a 
non- trivial motion, as predicted. 





(a) 



(b) 



Figure 11: A plane framework with Z 2 x Ci symmetry has a non-trivial flex on 
the fixed lattice. 



Here is the orbit matrix of the framework depicted in Figure 11 on the fixed 
lattice, with joints y>\, pi, Pz- 



{1,2} 
{2,3} 
U,3; 
{2,3; g 2 ] 
{1,3; g 3 ] 



Pi P2 P3 

(Pi ~ Pi) (P2 " Pi) 

(P2 " Pa) (P3 " Pi) 

(Pi " gl(Ps)) (p 3 " 

(p2 ~ gl(P3)) (P3 " g^iPlj) 

XP^-g^Pl)) (P2-^3 1 (P1)) 



In Table § we summarize the (Z 2 x C2)-symmetric Maxwell type counts 
for each of the lattice variants. For simplicity at this stage, we again assume 
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that no joint and no bar is fixed by the half-turn, so that all vertex orbits and 
edge orbits of the periodic orbit graph under the action of the group have the 
same size kc 2 = 2. For each type of lattice deformation, we always assume 
that e is chosen to be the least number of edges for the framework to be rigid 
without symmetry and to be compatible with the symmetry constraints given 
by Z 2 x C 2 . The number fe 2 in the final column denotes the dimension of the 
space of (Z 2 x C2)-symmetric infinitesimal flexes in each case. For 'generic' 
configurations, these extend to finite symmetry-preserving flexes. 



Table 3: Plane lattice deformations with C 2 symmetry. 



LatticeDef 
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ks 


ts 


4s 


rows 


columns -ts 


fs 


flexible 


c 2 


2 
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e = 2v + 1 


2v +3 


2 


distortional 


c 2 


2 





2 


e = 2v 


2v +2 


2 


hydrostatic 


c 2 


2 





1 


e =2v -l 


2v + l 


2 


fixed 


c 2 


2 








e =2v -l 


2v 


1 



6.2 Z 2 xC s - mirror symmetry in the plane lattice 

The mirror parallel to a side of the lattice restricts the possible lattices to rect- 
angles. This mirror symmetry is preserved by translation along the line of the 
mirror, so tc s = 1. 

We will consider periodic plane frameworks with symmetry Z 2 x C s again 
in two layers: (1) a fully flexible lattice; (2) a fixed lattice. 

Example 6.2.1: fully flexible lattice Z 2 x C s . The original (non-symmetric) 
necessary count for any minimally rigid periodic framework on the fully flexible 



lattice is e = 2v + 1 (recall Theorem 4.3 1. To permit mirror symmetry, with no 



vertex or edge fixed by the mirror, we will need to start with the shifted count 
2eo = 2(2i?o ) + 2 or equivalently eo = 2v§ + 1. 

Under the mirror symmetry with a flexible lattice which preserves the sym- 
metry, the orbit matrix has 2 columns under each orbit of vertices, plus £q, - 2 
lattice scaling columns for the mirror preserving flexes of the lattice. Since 
tc e = 1, we have the necessary symmetric Maxwell condition 

e > 2z; + 2 - 1 = 2v Q + 1 

for periodic rigidity. This inequality, together with the previous condition for 
minimal rigidity without the mirror symmetry, suggests that there is no added 



flexibility from this mirror symmetry. The example in Figure 12 illustrates such 
a situation with v$ = 3, eo = 7, and eo = 7 = 2co + 1- It is indeed rigid on a 
flexible lattice, up to vertical translation along the mirror line. □ 

Example 6.2.2: fixed lattice Z 2 x C s . The original necessary count for a periodic 
framework on the fixed lattice to be minimally rigid is e = 2v - 2. With added 
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gi = ((o,i), «9 

gi = ((0,0), a) 
g 3 = ((0,0), a) 

g 4 = ((-l,0),ff) 



(a) 



(b) 



Figure 12: The mirrors (vertical lines in (a)) fit only with the two scalings and 
this framework prevents those scalings. The orbit graph corresponding to this 
framework is shown in (b). 

mirror symmetry, we have e = 2eo and v = 2vq, so a minimally rigid orbit graph, 
realized with C s symmetry, will have 2eo = 2(2co) - 2, or eo = 2v$ _ 1- 

Under the mirror with a fixed lattice, the orbit matrix has 2 columns under 
each orbit of vertices. Moreover, we have Iq, = 1 since the translation along the 
axis preserves the mirror symmetry along with the periodic lattice. This creates 
the necessary symmetric Maxwell condition 



for periodic rigidity. Together with the count for minimal rigidity without 
symmetry, this suggests that there is no added flexibility from this symmetry. 
□ 

It turns out that for mirror symmetry, all of the variants of lattice deforma- 
tions produce no added motions. 

6.3 Table of groups for the fully flexible lattice in 2-dimensions 

Examples 6.1.1 and 6.1.2 indicate a process that can be applied to other plane 
symmetries which preserve the lattice. Each row in Table [I] presents the calcu- 
lation for a given plane wall-paper group which is presented as Z 2 » S. Recall 
that we are not including the plane wall-paper groups that have core glide 
reflections or 3-fold and 6-fold rotations, since they require some significant 
modifications of the simple pattern presented here (see also Section |&2) . Thus, 
we do not have 17 lines in the table. 



e >2v -l 
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In each row of Table[4j the calculation has several parts - each producing an 
integer: 

1. the number of edge orbits, eo, so that k^eo > 2{ksVo) + 1, which guarantees 
that we have at least the number of edges needed for a rigid periodic 
framework without symmetry. This means we need to add a modified 
constant For TableQ this value is always 1, and the number of rows 
is always eo = 2fo + 1 ■ 

2. t$ which is the dimension of the space of translations contained in the 
symmetry element of S. This will be 2 for the identity group, 1 for a single 
mirror, and otherwise. 

3. is which is the dimension of the space of lattice deformations which pre- 
serve all the symmetries in S or equivalently, the number of independent 
parameters in the lattice system (edge lengths and angles). 

4. the comparison of these numbers as the number of rows eo compared to 
the number of columns minus tg. 2vq + £$ - ts- 

5. the difference fs = 2v + €$ — ts - {2v + 1) = t s - t s - 1 which is the 
dimension of the guaranteed extra non-trivial motions of the symmetric 
framework, over the rigidity which the original count without symmetry 
promised. 



Table 4: The added flexibility induced by basic symmetries on a fully flexible 
2-D lattice for 1? x S. 
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2 
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// 
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*2222 


4 





2 
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1 
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c 4 


4 
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1 
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Civ 


4m 


*442 


4 





1 


2z> + l 
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6.4 Table of groups for the fixed lattice in 2-dimensions 

As we mentioned earlier, it can be of interest to consider periodic frameworks 
where the lattice is fixed. In the following table, each row will present the cor- 
responding calculation for a given plane wall-paper group which is presented 
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as Z 2 xS. As above, this analysis does not include rows for the hexagonal tiling 
or groups which include glide reflections. 

As in the previous section, in each row of Table|5j the calculation has several 
parts, each producing an integer. 

Table 5: The added flexibility induced by basic symmetries on a fixed 2-D lattice 
forZ 2 x<S. 
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ts 


rows 
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fs 
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Note that the number of rows (edge orbits), eo, is now such that kseo > 
2(ksVo) - 2. This guarantees that we have at least the number of edges needed 
for a rigid periodic framework on the fixed lattice without symmetry. This 
means we need to subtract a modified constant c = LjjjJ- For Table |5j c is 2, 1, 
or 0. 

Since for a fixed lattice, we clearly have = for each group S, the 
corresponding column is omitted in Table [5] 

Analogously to Table |4j the final column of Table [5] shows the difference 
fs = 2vq —ts~ (2vq — c) = c — ts which is the dimension of the guaranteed extra 
non-trivial motions of the symmetric framework, over the rigidity which the 
original count without symmetry promised. 

Note that Table [5] does not include all the point groups from Table |4] The 
groups we omitted only produce 0's in the last column. 

7 3-D periodic frameworks with symmetry: Z 3 x S 

We now apply the basic patterns of the previous sections to investigate the 
types of counts which arise for periodic structures with added symmetry in 
3-space. As happened in the plane, these symmetries can have three impacts: 

(a) the symmetry can restrict the possible shapes of the lattice cell or equiva- 
lently, the symmetry constraints leave a specific subset of £s flexes of the 
lattice structure which preserve the desired symmetry. 

(b) the symmetry can block some, or all, of the translations of the lattice 
structure, altering the basic count of ts; 

(c) the symmetry determines the order of the group, that is, the size ks of the 
orbits. 
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7.1 Z 3 x d - inversive symmetry in space 



Consider the inversive symmetry in 3-space with the center of symmetry at the 
origin. This operation (which in the Schoenfiies notation is called i) takes a joint 
p to a joint -p. In many tables of crystal symmetry this symmetry operation 
is called central symmetry and the crystals are called centrosymmetric. All 
shapes of lattices are possible, and these fit into the triclinic lattice system (three 
angle choices). In the Schoenfiies notation, if inversion is the only non-trivial 
symmetry operation, the group is written as C,. In the Hermann-Mauguin 
notation, it is written as 1, and in the orbifold notation, it is written as lx. 

As in the plane, if we have a center of inversion c, and a translation vec- 
tor t then there is another inversion centered at c + It. So, given the lattice 
of translations Z 3 and one center of inversion at the origin, there is a full lat- 
tice of inversions, with translations 5Z 3 , and the group of operations on the 



framework is written Z 3 » d (see also Figure 13 a)). 



Example 7.1.1: fully flexible lattice Z 3 x C,. The necessary count for any 
minimally rigid non-symmetric periodic framework on the fully flexible lattice 
ise = 3p+3 (recall Theorem |4.3| . To permit inversion symmetry (kc i =2) we will 
need to start with the shifted count 2eo = 3(2co) + 4 or equivalently eo = 3z>o + 2. 

Since the full flexibility of the lattice fits with the inversive symmetry, we 
still have lc, — 6- Further, when we move to the symmetric periodic orbit matrix 
under inversive symmetry, all of the infinitesimal translations disappear from 
the kernel, so that fp, = 0- This gives rise to the symmetric Maxwell condition 

eo ^ 3c + 6 

for periodic rigidity. The gap eo = 3c + 2 < 3z?o + 6 implies that a graph 
which counted to be minimally rigid without symmetry, realized generically 
with inversive symmetry on a fully flexible lattice now has a space of (finite) 
flexes of dimension 4. □ 

Example 7.1.2: fixed lattice 1? x C,. The necessary count for any minimally 
rigid non-symmetric periodic framework on the fixed lattice to be minimally 
rigid is e = 3c - 3. To permit inversive symmetry, we will need to start with the 
shifted count 2eo = 3(2vq) - lor equivalently eo = 3co - 1. 

Since we again have f^. = 0, we obtain the necessary symmetric Maxwell 
condition 

eo > 3v Q 

for periodic rigidity. The gap eo = 3co - 1 < 3co predicts a non-trivial finite 
flex in generic realizations with inversive symmetry on the fixed lattice. □ 
As a summary, here is the impact of inversive symmetry for each of the 



variants of lattice flexibility introduced in Section 4.3 
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(d) (e) 



Figure 13: In 3-D, one center of inversion repeats with half the period (a). An 
orbit framework with 2 orbits of vertices is shown in (b), with the group ele- 
ments associated with the directed edges listed in (c). Parts (d) and (e) illustrate 
building up the corresponding symmetric-periodic framework, moving from 2 
to 8 orbits of edges (d). 



Table 6: 3-D lattice deformations with C, symmetry. 
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7.2 Z 3 x Ci and Z 3 x C s - half -turn and mirror symmetry in space 



Assume we have a 2-fold rotational axis along the z direction. This places 
the pattern into the monoclinic crystal system: one face of the lattice is a 
parallelogram (perpendicular to the axis) and two faces are parallel to the axis 
and perpendicular to the parallelogram face. For this type of lattice, there are 4 
lattice parameters: the scale of each of the generating translations, and the one 
angle between the two generating translations of the parallelogram. 

Example 7.2.1: fully flexible lattice 1? x C 2 . With a fully flexible lattice, the 
necessary minimal number of edges for a periodic framework to be rigid and to 
be compatible with half -turn symmetry is 2eo = 3(2co) + 4, or eo = 3co + 2. In the 
orbit matrix, there are four columns corresponding to the lattice deformations, 
so the necessary symmetric Maxwell type count for periodic rigidity is 

e > 3y + 4 - 1 = 3z; + 3. 

Since we started with eo = 3vq + 2 < 3vq + 3, we predict a non-trivial symmetry 
preserving finite flex for generic realizations with half-turn symmetry on the 
flexible lattice. □ 

Example 7.2.2: fixed lattice Z 3 x C 2 . With a fixed lattice, the necessary minimal 
number of edges for a periodic framework to be rigid and to be compatible with 
half -turn symmetry is 2eo > 3(2co) — 2, or eo > 3vq - 1. The necessary symmetric 
Maxwell type count for periodic rigidity on the fixed lattice is 

e > 3v - 1. 

Thus, we do not detect any added motions in this case. □ 
In Table|7]we present the (Z 3 x C2)-symmetric Maxwell type counts for each 
type of lattice deformation. 



Table 7: 3-D lattice deformations with C 2 symmetry. 
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Consider a periodic framework in space with mirror symmetry. For this 
new group, there are only two key calculations to be done: 

1. tc s = 2, since the two translations on directions within the mirror will 
(instantaneously) preserve the mirror. 
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2. €c s = 4. Although there initially appear to be two alignments for the 
mirror: (i) parallel to two translation axes and perpendicular to another 
or (ii) containing an axis of translation, these turn out to be two variations 
of the same larger space tiling, and crystallographers only consider the 
first version. In this case we have an orthorhombic lattice system, and we 
have four parameters, (q s — 4. 



Example 7.2.3: fully flexible lattice Z 3 x C s . As before, we start with the 
following initial count without symmetry: 2eo = 3(2vq) + 4, or eo = 3co + 2. 
From the periodic symmetric orbit matrix we obtain the following necessary 
symmetric Maxwell type count for periodic rigidity: 

e > 3v Q + 4 - 2 = 3v + 2. 

This suggests that there is no additional flexibility in the structure when mirror 
symmetry is added. □ 
It turns out that for mirror symmetry, all of the variants of lattice deforma- 
tions produce no added motions. 



7.3 Table of groups for the fully flexible lattice in 3-dimensions 

Following the process illustrated in the previous examples, we can track the 
necessary increases in flexibility which follow from minimal generically rigid 
periodic frameworks for various symmetry groups Z 3 »S in 3-space. As before, 
this does not include rows for the groups with 6-fold rotational symmetry, or 
any patterns with glide reflections. They will require some significant modifi- 
cations of the simple pattern presented here. 

Analogous to the tables in Sections 6.3 and 6.4 in each row of Table |8j the 



calculation has several parts - each producing an integer. 

The number of rows (edge orbits), eo, is such that k$eo > 3(ksz>o) + 3, which 
guarantees that we have at least the number of edges needed for a rigid periodic 
framework without symmetry. This means we need to add a modified constant 
c = [~^~|. For Table|8j c = 3 for ks = 1, c = 2 for k$ = 2, and c = 1 for all bigger 
orbit sizes. 

As usual, t$ is the dimension of the space of translations contained in the 
symmetry element of S. This will be t$ — 3 for the identity group, t$ = 2 for 
a single mirror, t$ = 1 for a single rotation (with or without a mirror along the 
axis), and tg = if only a point is fixed. 

In Table|8]we compare the number of rows, eo, with the number of columns 
minus t s , 3v Q +£s~ts} the difference fs = 3c + £s _ ts ~ (3^o + c) = (s ~ is ~ c 
is the dimension of the guaranteed extra non-trivial motions of the symmet- 
ric framework over the rigidity which the original count without symmetry 
promised. 
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Table 8: The added flexibility induced by within basic symmetries on a fully 
flexible 3-D lattice for 1?»S. 
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7.4 Table of groups for the fixed lattice in 3-dimensions 

In Table [9j we track the necessary increases in flexibility which follow from 
minimal generically rigid periodic frameworks on a fixed lattice for various 
symmetry groups in 3-space. This analysis is analogous to the one in the 
previous section. We simply remove the column for ig which is always 0, and 
work with the modified counts. 

The entries (-1) in Table [9] indicate that, for this group, the symmetry guar- 
antees that there is a symmetric self-stress in the symmetric framework (see 



also Section 8.5 1. Because the patterns of and occasional (-1) become clear 



quickly, we do not fill in all rows of the matrix. 
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Table 9: The added flexibility induced by symmetries on a fixed 3-D lattice for 



Lai. System 


Schs 


H-M s 


orbs 


ks 


*5 


rows 


columns -ts 


fs 


triclinic 


C x 


1 


11 


1 


3 


357-3 


3v-3 







Ci 


1 


lx 


2 





3v - 1 


3v -0 


1 


monoclinic 


c 2 


2 


22 


2 


1 


3^o - 1 


3v - 1 





// 


C s 


m 


1* 


2 


2 


3v - 1 


3v -2 


(-1) 


// 


C211 


2/m 


2* 


4 





3v 


3v -0 





ortho rhomb. 


Civ 


222 


*22 


4 


1 


3v Q 


3v - 1 


(-1) 


// 


£>2 


mm2 


222 


4 





3v Q 


3v -0 





// 




mmm 


*222 


8 





3v q 


3v -0 





tetragonal 


C 4 


4 


44 


4 


1 


3v 


3v - 1 


(-1) 




S 4 


2 


2x 


4 





3v Q 


3v -0 






8 Extensions and further observations 
8.1 Adjusting for fixed joints and bars 

To keep the analysis simpler, up to this point, we have not considered joints or 
bars which are fixed by a non-trivial symmetry operation in S. This was for 
simplicity of our counts (keeping a fixed ks for all vertex and edge orbits), but 
incorporating these modifications are not a barrier to the final analysis. The 
methods presented here also apply to symmetric periodic frameworks with 
fixed structural components, as has been shown in the previous work on finite 
symmetric frameworks ( jSchulze and W hiteley |2010) . We give a brief indication 



of how this can work, in two cases. 

We first consider the case where the given periodic framework has a bar 
(but no joint) which is fixed by a non-trivial symmetry operation in the group 
S. In this case, the comparison of counts needs to be adjusted, since there are 
now edge orbits of different sizes - the size of an edge orbit corresponding to a 
bar which is fixed by a non-trivial symmetry operation will no longer be equal 
to the order of the group S. This is illustrated in the following example. 

Example 8.1.1 The necessary count for any minimally rigid non-symmetric 
periodic framework on the fully flexible lattice is e = 3v + 3. If a framework 
with Z 3 x d symmetry has exactly one bar which is fixed by the inversion, 
then we have e = 2e§ - 1 . Thus, in this case, we do not need to shift the count 
e = 3v + 3 to allow for inversive symmetry. We have 2eo — 1 = 3(2co) + 3, and 
hence eo = 3c + 2. The symmetric Maxwell condition for periodic rigidity is 
eo ^ 3^o + 6. So we detect 4 added degrees of flexibility. □ 
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The method for adjusting all of the table entries, for any number of fixed 
bars, should now be accessible to the reader. 

Similarly, for a periodic framework with symmetry Z d x S, there can also be 
joints in the unit cell which are fixed by some non-trivial symmetry operations 
in the group S. For these, the methods applied in Schulze and Whiteley (2010) 



also immediately transfer. Note, however, that if there exist joints that are fixed 
by non-trivial symmetry operations in S, we may not only have vertex orbits 
of different sizes, but the sets of columns corresponding to the vertices in the 
orbit matrix may now also be of varying size. 

Depending on what subspace the joint is now restricted to, the number 
of corresponding columns in the orbit matrix is reduced. If a joint of a 3- 
dimensional structure, for example, is restricted to a mirror plane, the number 
of columns will be reduced to 2. If it is restricted to a line (e.g. the intersection 
of two mirrors, or an axis of rotation) there will only be 1 column, and if it is 
restricted to a point (e.g. the center of inversion, or the intersection of a mirror 
and an axis) the number of columns will be 0. 

We give two samples to illustrate this. In both cases we assume that there 
are no bars which are fixed by a non-trivial symmetry operation in S. 

Example 8.1.2 The necessary count for any minimally rigid non-symmetric 
periodic framework on the fully flexible lattice is e = 3v + 3. If a framework 
with Z 3 x C s symmetry has exactly one joint which is fixed by the reflection in 
C s , then we have v = 2vq - 1 . Thus, in this case, we do not need to shift the count 
e = 3v + 3 to allow for mirror symmetry. We have e = leg - 3(2vq - 1) + 3 = 6vq, 
and hence - 3vq. The symmetric Maxwell condition for periodic rigidity is 

e > 3(v - 1) + 2 + 4 - 2 = 3v + 1. 

So we detect one added degree of flexibility. □ 

Example 8.1.3 As noted above, the necessary count for any minimally rigid 
non-symmetric periodic framework on the fully flexible lattice is e = 3v + 3. 
Recall from Example 7.2.3 that we need to shift this count to e = 3v + 4 if we 
want to permit Z 3 x C s symmetry, where no joint and no bar is fixed by the 
reflection in C s . No added motion was detected in this case. We now start with 
a framework whose underlying periodic graph satisfies the count e = 3v + 4, 
and we further assume that exactly two of the joints are fixed by the reflection 
in C s (and hence lie on two distinct points on the corresponding mirror plane). 
Then we have v = 2vo - 2. Thus, we have e = 2eo = 3(2vo - 2) + 4 = 6vq - 2, and 
hence eo = 3^0 ~ 1- The symmetric Maxwell condition for periodic rigidity is 

eo>3(z> -2) + 2(2)+4-2 = 3z;o- 

So we now detect one added degree of flexibility, although the framework is 
overbraced, generically. □ 
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8.2 Other groups 



Given the wide array of lattice systems and space groups which occur in crys- 
tals, we have only analyzed some of the possibilities. The groups we considered 
were restricted to the form Z d x S. 

It is a small change to include groups with glide reflections as genera- 
tors which are not semi-direct products. Such a group does include Z d as a 
subgroup, so it admits a representation as a periodic orbit graph. The glide re- 
flection acts within the periodic orbit graph, and the same counts and methods 
presented here can be applied (see also Example 8.1.3). Note, however, that the 
group acting on the periodic orbit framework is not a point group (i.e., a group 
of isometries leaving a point fixed), so it is technically distinct from the groups 
we have considered so far. 

Example 8.2.1 We start with the necessary count for minimal rigidity of a non- 
symmetric periodic framework on the fully flexible lattice: e = 3v + 3. This 
needs to be shifted to 2e = 3(2v ) + 4 or e = 3v + 2 to permit glide reflectional 
symmetry. As in the case of Z 3 x C s , there are 4 parameters for the lattice 
deformations in the orbit matrix, and there are still two translations which 
preserve the glide reflection, along with the periodic lattice. This leads to the 
following symmetric Maxwell condition for periodic rigidity: eo > 3vo + 4 - 2 = 
3v + 2. So, as in the case of Z 3 x C s , we do not detect any flexes with these 
counts. 

Note that placing a joint on the mirror plane of the glide reflection does 
not mean that this joint is fixed by the glide reflection, but it does restrict the 
number of degrees of freedom for this joint - that is, the number of columns 
corresponding to this joint in the orbit matrix - to 2. Thus, shifting two joints 
onto the glide plane will generate an added degree of flexibility, as it did for a 
single mirror (recall Example 8.1.3). □ 

It is also simple to include a single screw symmetry in 3-space, such as a 2- 
fold screw. This will have the same impact as the 2-fold rotation, generating an 
additional flexibility in the flexible lattice. However, we have not completed the 
analysis for more general mixtures of glide reflections, and screw symmetries, 
with the various compatible lattices and groups. We anticipate that it will 
be possible to extend all of the results to the space groups which share the 
underling lattice systems we explored here. 

We have also not analyzed the symmetries which include 6-fold rotations 
in the plane or in 3-space - those related to the triangular lattice. These triangle 
groups require some additional care with the gain graphs and the orbit counts 
when they are placed into the unit cells as we factor out the translations. Again, 
we anticipate these lines of the tables will fall into place once this analysis is 
completed. In short, we anticipate that it will be possible to give analogous 
orbit matrices and counts for all of the plane wallpaper groups and all of the 
230 space groups. 

However, the class already covered includes a number of key examples, 
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and these example do illustrate that symmetry can add flexibility to structures 
which otherwise count to be generically rigid. 



8.3 Higher dimensions 

The basic results for rigidity matrices for periodic structures and for orbit 
rigidity matrices for symmetric structures all extend to arbitrary dimensions 
( |Borcea and Streinu[|20~10b|a{|Schulze and Whiteley[|2010} . 

All of the techniques for combined analyses also extend to arbitrary dimen- 
sions d for groups of the form Z d x S. As an example, we summarize the 
extension for the groups Z d x C, and Z d x C s . 

For (non-symmetric) generic rigidity on a full y fle xible lattice in d-space, a 
necessary condition is e > dv + (^) (recall Theorem 4.3 1. In both of the following 
examples, we will assume that we have the minimal number of edges for the 
graph to be rigid on a fully flexible lattice and to be compatible with the given 
symmetry, with no joint and no bar fixed by a non-trivial symmetry operation. 

Example 8.3.1 Consider a periodic framework with inversive symmetry in 
dimension d and the associated orbit matrix. We have 2eo = d(2vg) + (^) if (2) 
is even, and 2eo = d(2vo) + u) + 1 otherwise. Further, we have the following 
column count for the orbit matrix: dz>o + ( t 1 ). There are clearly no residual 
translations which preserve the inversive symmetry along with the periodic 
lattice. So, since 



e = dv + 



/2 



< dvo + 



d + 1 



we detect fc t = f^ 1 ) _ [(2)/^] degrees of flexibility. This is a number which is 
growing as [.(^/^J + d, which is quadratic. □ 

Example 8.3.2 For mirror symmetry, we also have a pattern across the dimen- 
sions. In dimension d, there will be one lattice vector perpendicular to the 
mirror, and the rest will be parallel to the mirror. This removes a set of d — 1 
angles as lattice variables, leaving £q s — m + 1. There are also d — 1 residual 
translations. So, since 



e = dv + 



< dv + 1 I + 1 - (d - 1), 



we detect f Cj = (^) + 1 - (d- 1) - [(2V2] degrees of flexibility. This simplifies as 
fc, - L(^) /2J + 2 - d. For d > 3 this is positive and growing quadratically. □ 
Similar methods can be used to analyze a number of other groups of the 
form Z d x S in higher dimensions, provided we have worked out the flexibility 
of the lattices which support the symmetry, in that space. 
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8.4 Topology of orbit frameworks 



We may regard orbit frameworks as graphs embedded on appropriately cho- 
sen orbifolds (generalizations of manifolds, which locally resemble Euclidean 
space). This orbifold is defined by the original setting of the framework (R 2 
or M 3 ), modulo the symmetry group. For example, periodic frameworks have 
symmetry group Z d , and may be viewed as graphs on the d-dimensional topo- 
logical torus R d /Z rf . We can think of this as the unit cell with pairs of opposite 
faces identified. Similarly, a plane framework with C n symmetry can be re- 
garded as a framework on a cone, with opening angle 2n/n. The orbit matrix 
now also provides conditions for rigidity and flexibility for the (orbit) frame- 
works on this surface. 

For periodic frameworks with additional symmetry, the underlying orbifold 
may be more exotic. For example, periodic frameworks with mirror symmetry 
in the plane or space (given by the groups Z 2 x Cs and Z 3 x Cs respectively) 
correspond to frameworks on 2- or 3-spheres S 2 and S 3 , but with a flat metric. 
Frameworks in 3-space with inversive symmetry (Z 3 x C,) have an orbifold 
with topology of P 3 , projective 3-space. Similar statements are possible for 
all frameworks which admit an orbit framework under the action of their 
symmetry group. Again, the periodic symmetric orbit matrices represent the 
rigidity matrices for frameworks actually living in these more exotic spaces, 
with flat metrics. The results here give some necessary conditions for rigidity 
on these orbifold surfaces. 

As an additional topological and geometric layer, the more detailed stud- 
ies of the orbifolds for space groups describe a number of these orbifolds as 
fiberfolds |Conway et al. 2001} . These are essentially 'fibered prisms' over the 
plane orbifold. For a group S in the plane, we have the stretched group S such 
that s(x, y, z) = (s(x, y), z) for each group element. The simplest extension of 
our counts occurs for the scaling lattice deformations. In this setting, we have 

= ks, = ts + 1 and €g=€s + l. Comparing eo = 2vq for the plane without 
symmetry with the scaling counts eo = 2co + £g — t$, we have a flexibility of 
is - is- Doing the same calculation in 3-space, we find - tg = is ~ is- We 
conclude that the flexibility induced by S and by S is the same. 

In general, it is an open problem to investigate this connection across di- 
mensions more thoroughly. The examples in 6 8.3 give a hint of how some of 
this might work. It is also an interesting problem to predict how flexibility of 
a periodic symmetric graph in one dimension connects to flexibility for some 
associated periodic symmetric graph in the next dimension up. 



8.5 Statics and stresses 



In this paper, we have not focused on stresses - row dependencies of the (pe- 
riodic) rigidity matrices, or symmetric stresses - row dependencies of the orbit 



matrices (Schulze and Whiteley 2010 White and Whiteley 1983 Whiteley 
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19961. Imposing symmetry on a periodic structure, however, may not only 



increase the flexibility of the structure, but it can also increase the dimension of 
the space of self -stresses of the structure (recall Table |9j for example). There are 
essentially three possible sources for these additional self-stresses. 

As a first source of stress, in a number of cases, we had to have additional 
edges to ensure the count of edges was divisible by kg. This of course leads 
to additional row dependencies of the (periodic) rigidity matrix which can be 
detected by simply counting the number of vertices and edges of the under- 
lying periodic orbit graph of the structure. These stresses may or may not be 
symmetric. 

Secondly, where introducing symmetry takes a minimally rigid periodic 
framework to an infinitesimally flexible framework this is a guarantee that the 
symmetric periodic framework has a space of self -stresses of dimension equal 
to the dimension of the added infinitesimal flexes. Note that if we detect a 
symmetric infinitesimal flex in the periodic framework from the orbit counts, 
then the framework must have an asymmetric self -stress. Similarly, if we detect 
a symmetric self-stress in the periodic framework from the orbit counts, then 
the framework must have an asymmetric infinitesimal motion. This is because 
the orbit counts cannot detect the presence of 'paired' symmetric infinitesimal 
flexes and self stresses. 

The third possibility is that although our orbit counts do not detect any 
self-stresses, the rows of the orbit matrix are not actually independent. To 
determine this requires additional direct computation of the rank of the orbit 
matrix. 

The classic plane example of Bottema's mechanism - K44 with symmetry 
group C20 in the plane - is known to be flexible not from the simple orbit counts, 
but from these counts plus the added information that the orbit matrix has a row 



dependency (Schulze and Whiteley 20101. One can make periodic structures 



from this example, and create added flexes due to the periodic symmetry, which 
are not predicted simply by the counts presented here for the rows and columns 
of the orbit matrix. 



8.6 Sufficient conditions for rigidity or independence with 
symmetry 

In the classical work for finite frameworks, there are results where combinatorial 
properties on a graph, and its subgraphs, guarantee that the rigidity matrix is 
full rank, or independent, at generic configurations. The most famous example 
in the plane is Laman's theorem (Laman |1970 Whiteley 1996) . In 3-space, 
some partial sufficiency results come from graphs constructed by specific types 
of inductive constructions, or planar graphs which are 3-connected < |Whiteley] 
1996) . 

One can ask about combinatorial properties which guarantee that a periodic 
orbit matrix has full rank. This has been determined in the plane for both 
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the flexible lattice (Malestein and Theran 20101 , and the fixed lattice ( |Ross 



2010| . In addition, for some plane symmetry groups, there are combinatorial 



characterizations of graphs which are minimal rigid and maximal independent 
frameworks at symmetric generic configurations, as well as conjectures for 



other cases < |Schulze 2010c) . It is then a natural question to be explored to find 



combinatorial characterizations of graphs which are minimal periodic rigid 
with a given symmetry, in the plane. This may be accessible for at least some 
cases. 

One can also ask for characterizations of special classes of symmetric pe- 
riodic rigid graphs in 3-D, for example, frameworks built by connecting sym- 
metric convex polyhedra in a periodic form, or frameworks built by inductive 
constructions based on edge and vertex splitting. 

In summary, our initial exploration presented here raises a large number of 
interesting unsolved problems, and areas for further research. We invite the 
reader to join in this exploration. 
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